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ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ТЕМПЕРАТУРНЫХ ПОЛЕЙ СИЛОВЫХ ТРАНЗИСТОРОВ 

С УЧЕТОМ РАЗРЫВОВ КОЭФФИЦИЕНТОВ ПЕРЕНОСА 
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Численно решена двумерная задача о нестационарном температурном поле в корпусе силового биполярного транзистора с во- 
семью источниками тепловыделения общей мощностью 60 Вт. Исследовано влияние сеточных параметров на точность вычисле- 
ния температур. Установлено, что при определенных сеточных параметрах возможно достижение очень высокой точности вычис- 
лений. 


Введение 

Выбор условий работы современных полупро- 
водниковых приборов (ППП) в том числе транзис- 
торов с большим энерговыделением во многом оп- 
ределяется тепловыми режимами их эксплуатации. 
Экспериментальная отработка оптимальных (хотя 
бы по отдельным критериям) тепловых режимов 
представляет очень сложную задачу в связи с малы- 
ми размерами ППП, малыми временами работы, 
высокими значениями рассеиваемой мощности и 
значительными скоростями роста и падения темпе- 
ратур элементов ППП. Измерения температур на 
отдельных поверхностях таких элементов возмож- 
ны только со значительными погрешностями из-за 
высокой инерционности как термопарных, так и 
оптических методов измерения. Оценка же уровня 
температур на закрытых границах между элемента- 
ми (или внутри какого-либо элемента) практически 
невозможна потому, что ввод термопар в зону изме- 
рения существенно трансформирует поле темпера- 
тур в этой зоне и нарушает условия работы ППП. 

В связи с этим оценки диапазонов изменения 
температур внутри полупроводникового прибора с 
приемлемой для проведения проектных и конструк- 


торских работ точностью возможны только по ре- 
зультатам математического моделирования процес- 
сов переноса теплоты с учетом реальных механиз- 
мов теплообмена на всех границах [1,2]. Известны 
результаты такого моделирования для приборов с 
относительно низкой рассеиваемой мощностью 
(ед. Вт) и достаточно большими характерными вре- 
менами работы (сотни с) [3, 4]. Для ППП типа мощ- 
ных транзисторов, работающих менее секунды, ка- 
ких либо результатов численного моделирования до 
настоящего времени не опубликовано. 

Целью данной работы является решение задачи 
о нестационарном температурном поле мощного 
транзистора в рамках двумерной модели теплопро- 
водности для области с разрывными значениями 
теплопроводности, теплоемкости и плотности. 

Постановка задачи 

Анализ типичных конструктивно компоновоч- 
ных решений по схемам мощных транзисторов пока- 
зывает, что характерными для таких приборов явля- 
ются сильные разрывы (в 100 и более раз) коэффици- 
ента теплопроводности на границах между конструк- 
тивными элементами. Значения теплоемкости и 
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плотности на границах также изменяются, но не бо- 
лее чем 6 раз. Эти особенности рассматриваемого ти- 
па задач ранее специально не исследовались. Поэто- 
му нет отработанных алгоритмов их решения. 


Рассматривался типичный вариант исполнения 
полупроводниковых приборов в корпусе ТО-220. 
За основу взят широко распространенный транзис- 
тор КТ819Г. Транзистор выполнен в металлоплас- 
тиковом корпусе рис. 1, 2. Основой корпуса являет- 
ся медное основание (зона 1, 3) толщиной 1,22 мм, 
на котором припаян полупроводниковый кристалл 
(зона 4), пластик закрывает кристалл и часть мед- 
ного основания (зона 2,3,4). 




Рис. 1. Общий вид транзистора КТ819Г в корпусе ТО -220: 
а) фотография; б) схема 


В кристалле находятся 8 транзисторов. Элект- 
роды соединены параллельно на поверхности 
кристалла медным напылением толщиной 15 мкм. 
Транзистор имеет тепловой контакт с радиатором. 



Рис. 2. Схема области решения; 1) медное основание, 
2) пластик, 3) пластик на медном основании, 
4) кристалл на медном основании в пластике, 5) ба- 
зовое медное напыление на кристалле, 6) эмиттер- 
ное медное напыление на кристалле, 7) заштрихо- 
ванные области восьми источников тепловыделения 


При постановке задачи приняты следующие до- 
пущения: 

1. В медном основании (зона 1) не учитывается 
наличие крепежного отверстия. Численный 


эксперимент показывает, что градиент темпера- 
туры в этой зоне мал, поэтому отсутствие отве- 
рстия в этой зоне слабо влияет на распределе- 
ние температур. 

2. Не учитываются выводы транзистора. Посколь- 
ку выводы вплавлены в пластик с низкой теп- 
лопроводностью, а к кристаллу припаяны тон- 
кими проводниками, выводы не имеют хороше- 
го теплового контакта с областью интенсивного 
тепловыделения. 

3. Интенсивность источника тепловыделения не 
зависит от координаты. 

4. Коэффициенты теплоотдачи и теплопередачи 
не зависят от температуры. 

5. Температура радиатора постоянна и равна тем- 
пературе внешней среды. 

Численный анализ температурных полей про- 
водился при теплофизических характеристиках, 
приведенных в табл. 1. 


Таблица 1. Теплофизические характеристики материалов 


Материал 

Теплоемкость, 

Дж/(кг*К) 

Плотность, 

кг/м 3 

Теплопроводность, 
Вт/(м -К) 

Медь 

390 

8930 

365 

Компаунд 

2000 

1500 

0,8 

Кремний 

712 

2300 

159 


Поскольку толщина элементов прибора зависит 
от координат хи у а задача решена в двумерной 
постановке, то все размеры по высоте приведены к 
одной величине. Вследствие этого изменятся экви- 
валентные теплофизические характеристики неко- 
торых областей. В данном случае все вертикальные 
размеры приведены к толщине медного основания 
4= 1,22 мм. С учетом такого усреднения значения 
теплофизических характеристик определяются из 
выражений (1-3): 

N 

2>,ч 

К=^—, ( 1 ) 

где А„ 4 - коэффициент теплопроводности и толщи- 
на /-ого слоя вертикальной структуры соответствен- 
но, А 0 - эквивалентный коэффициент теплопровод- 
ности всех слоев УѴ, приведенный к толщине 4- 

N 

X А • 4 

Ро=^Ч > ( 2 ) 

а о 

где Рі - плотность /-ого слоя вертикальной структу- 
ры, р 0 - эквивалентная плотность всех слоев УѴ, 
приведенная к толщине 4- 

N 

X С І ■ А ■ 7 

„ 2=1 /^\ 


где Сі -удельная теплоемкость /-ого слоя вертикаль- 
ной структуры, с 0 - эквивалентная удельная тепло- 
емкость всех УѴ слоев, приведенная к толщине 4- 
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Коэффициенты конвективной теплоотдачи на 
границах х=0; х=Ь х ; у=0; у=Ь у приняты равными 

ос 2 = 2 . Коэффициент теплопередачи через 

м -К 

поверхность контакта с радиатором, принят рав- 
ным а х = 80 кВт , что соответствует толщине теп- 
м -К 

лопроводной пасты 10 мкм и ее удельной теплоп- 
роводности 0,8 -Щт. 

м • К 

Математическая модель 

В такой постановке задача сводится к решению 
нестационарного уравнения теплопроводности (4) 
с соответствующими граничными (6-9) и началь- 
ными (5) условиями: 

ч , ч дТ д дТ 

С(х, у) ■ р(х, у)- — = — (Цх, у)-—) + 

ді ох ох 

д дТ а х • (Т ѵ - Т) / 4 ч 

+ т - (^( х > у) ' ~) + і 2( х > у) + 5 ^ 

ду ду а 0 

где Т - температура, С(х,у) - удельная теплоем- 
кость, р(х,у) - плотность, Цх,у) - теплопровод- 
ность, ()(х,у) - интенсивность тепловыделения, 
й - приведенная толщина прибора. 

0</<оо, 0<х<Ь х , 0 <у<Ь у , 

где Ь х , Ь у - максимальные линейные размеры. 


1 = 0; Т = Т 0 (х,у); 

(5) 

х = 0; -хЖ = а 2 -(Т ѵ -Т); 

(6) 

х = Т,\ Ь-%=<*ЛТ Ѵ Т); 

(7) 

У = 0; А • ^ = а 2<Т, Т); 

(8) 

у = Д; А • ^ = а 2 • (Т ѵ Т ). 

(9) 


Задача решена методом конечных разностей, 
разностные аналоги уравнения решены методом 
прогонки [2]. 

Результаты численного исследования 

Численный анализ выполнен при суммарной ин- 
тенсивности тепловыделения 60 Вт, равномерно расп- 
ределенной по всем восьми локальным источникам. 

Численное моделирование велось на неравно- 
мерной сетке. Шаг сетки по координатам хи у из- 
менялся от некоторого минимального значения к тп 
в геометрической прогрессии с коэффициентом 
1,3. Шаг сетки по координатам увеличивался от 
внешней границы элемента до середины области 
решения и затем уменьшался к противоположной 
границе элемента. Таким образом, сетка сгущалась 
на границах между элементами. Численные иссле- 
дования велись в диапазоне изменения времени до 
3 с. Начальная температура принималась равной 
температуре окружающей среды 300 К. Шаг по вре- 
мени составлял от 0,1 до 12,8 мкс. Минимальный 
шаг по координате к тт варьировался от 1 до 45 мкм. 

Типичное распределение температурного поля 
при /=3 с, минимальном шаге по координате 
к Ып = 32 мкм и шаге по времени т=2Т0 _6 с приведено 
на рис. 3. 

Хорошо видны максимумы температур, соотве- 
тствующие зонам интенсивного тепловыделения 
транзистора. Максимальные перепады температур 
достигают в отдельных сечениях 23 К. Соответ- 
ственно отклонение от средних значений составля- 
ет ±11,5 К. Анализ температурного поля, приве- 
денного на рис. 3, показывает, что оно в целом хо- 
рошо соответствует условиям тепловыделения и 
теплопереноса в рассматриваемой структуре. При 
этом необходимо учитывать, что рамки двумерной 
модели несколько ограничивают возможности 


т, к 
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Рис. 3. Распределение температур на 3-й секунде работы ППП ( транзистора) 
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численного анализа. Если ориентироваться на ре- 
зультаты [3, 4], то можно предположить, что пере- 
ход к пространственной модели теплопереноса в 
силовом транзисторе может позволить уточнить 
полученные результаты. Но характер установлен- 
ных закономерностей при этом не изменится. 


Решение же пространственной задачи теплопере- 
носа для рассматриваемого объекта исследования 
возможно только на сетке, параметры которой выбра- 
ны, если и не оптимальными, то по крайне мере наи- 
более рациональными и обоснованными из всего 
многообразия возможных комбинаций сеточных па- 
раметров. До настоящего времени теоретически обос- 
нованные рекомендации по практической реализа- 
ции такого выбора для областей с сильными разрыва- 
ми коэффициентов переноса отсутствуют. По край- 
ней мере, для исследуемых в данной работе объектов. 
Поэтому специально решалась задача выбора сеточ- 
ных параметров, обеспечивающих максимально воз- 
можную точность вычислений температур с учетом 
особенностей рассматриваемой модели теплоперено- 
са. В связи с отсутствием экспериментальных данных 
и невозможностью их получения в реальных услови- 
ях, оценка достоверности полученных в ходе вычис- 
лительных экспериментов результатов по температур- 
ным полям силового транзистора проводилась на ос- 
новании анализа точности выполнения закона сохра- 
нения энергии. При реализации каждого варианта 
численного решения задачи рассчитывалась погреш- 
ность баланса энергии в области решения. Также 
контролировались максимальные отклонения темпе- 
ратур и средние отклонения температур от поля тем- 
ператур варианта с самым мелким шагом. Погреш- 
ность по балансу энергии определялась по формуле 


8 = 


А-(Ѵ х+ Ѵ 2+ Ѵ ъ ) 


А 


100 %; 


4 = 


N Мх Му 

І2І«. 

п= 1 х=1 у = 1 


(С-Г)А.А.^, 


у- 


где Ѵ\ - тепло, отведенное от корпуса прибора че- 
рез нижнюю поверхность за время УѴт, М х и 
М у - число узлов координатной сетки по х и у соот- 
ветственно, 7Ѵ- число шагов по времени. 


N Мх 

77 = 1 Х = 1 


Г\ +Т п м -2-Т 

х,\ х,Му ѵ 


-к - т + 


где А - суммарное тепловыделение всех 8-ми ис- 
точников. 


АТ 


шах 


тах(7 ] 




где, х=1 ...М х , у=1...М г 



Мх Му 



х=1 у = 1 


Т -Т 

1 Х,у 2 х,у 





При варьировании шага по времени определя- 
лась погрешность по балансу энергии, максималь- 
ное отличие температуры по (10) и среднее отличие 
температуры по (11) от случая с самым мелким ша- 
гом по времени т=210 _6 с. В табл. 2 приведены ре- 
зультаты вычисления погрешности расчетов при 
различных шагах по времени при к тіп = 32 мкм, 
М= 136, 4^=156. Варианты сравнивались с базо- 
вым вариантом при по максимальным (10) и сред- 
ним (11) отличиям температур. 


Таблица 2. Зависимость погрешности вычисления баланса 
энергии и отклонения температур от шага по вре- 
мени 


т, мс 

Е, ІО' 3 % 

т т , К 

ІО" 3 К 

0,002 

0,048 

0 

0 

0,1 

2,386 

0,052 

2,054 

0,2 

4,781 

0,102 

4,495 

0,4 

9,574 

0,219 

9,035 

0,8 

19,17 

0,452 

18,115 

1,6 

38,37 

0,919 

36,275 

3,2 

76,82 

1.853 

72,585 

6,4 

153,9 

3,721 

144,72 

12,8 

308,8 

7,465 

290,52 


Хорошо видно, что уменьшение шага по време- 
ни ведет к снижению погрешности по балансу 
энергии в системе. Причем приемлемые результа- 
ты получаются только при значениях г в десятые 
доли мс. При этом отклонения А Г тах могут состав- 
лять всего менее 0,1 К, при погрешности 8 соотве- 
тственно 0,003. ..0,004 %. Полученные результаты 
показывают, что даже при разрывных теплофизи- 
ческих характеристиках с большими отклонения- 
ми Я на границах между элементами транзистора 
может быть достигнута приемлемая точность вы- 
числений за счет использования очень мелких ша- 


N Му 
>7=1 у = 1 


гг, 11 . грп 

м ,У ^ Мх,у 


-2-Т 


К ' т > 

у У 


где ІІ 2 - тепло, отведенное от корпуса прибора че- 
рез боковые поверхности за время ТѴ-т. 


V, = 


Мх Му 

22ч, 

х=1 у = 1 


- р - (Т -Т) - сі- И -к 

' Х,у V х,у V / Х х у 


где ІІ 3 - тепло, накопленное внутри прибора за счет 
теплоемкости. 


Мх Му 


А = 



х=1 у = 1 


-сі -к -к -N -т, 

Х,У * х У у ’ 


гов по времени. 

Таблица 3. Зависимость погрешности вычисления баланса 
энергии от минимального шага по пространствен- 
ной координате при т=0,4 мс 


4™,. мкм 

и 

ч 

Е, ІО' 3 % 

1 

460 

652 

9,622 

2 

380 

536 

9,628 

4 

316 

422 

9,614 

8 

252 

324 

9,613 

16 

188 

288 

9,597 

32 

136 

156 

9,574 

45 

116 

116 

9,577 
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Рис. 4. Распределение разности температур при изменении шага по времени от 0,2 до 0,4 мс 


На рис. 4 приведены разности температур Т 3 двух 
вариантов расчетов (при шаге по времени 0,4 и 
0,2 мс) по формуле 

т - Т - Т 

Ъх,у 2 х,у \х,у ' 

Видно, что эти разности не превышают 0,12 К. 
Но также видно, что максимальные отклонения 
достигаются в зонах наибольших тепловыделений, 
расположенных вблизи границ между элементами. 

В табл. 3 приведены результаты вычисления 
погрешности баланса энергии в зависимости от 
минимального шага сетки по координатам, при 
шаге по времени т=0,4 мс. При варьировании ми- 
нимального шага по координате погрешность оп- 
ределялась только по балансу энергии, поскольку в 
каждом случае узлы неравномерной сетки имели 
несовпадающие координаты. 


СПИСОК ЛИТЕРАТУРЫ 

1. Дульнев Г.Н. Тепло- и массообмен в радиоэлектронной аппа- 
ратуре. — М.: Высшая школа, 1984. —247 с. 

2. Самарский А.А. Теория разностных схем. — М.: Наука, 1977. — 
656 с. 

3. Алексеев В.П., Кузнецов Г.В., Шлома С.В. Расчет температур- 
ного поля печатной платы с учетом конвективного и радиаци- 


Выводы 

Результаты численного эксперимента показы- 
вают, что погрешность численного решения ли- 
нейно связана с шагом по времени и слабо зависит 
от шага по координате. Наибольшая погрешность 
возникает на границах неоднородностей с больши- 
ми перепадами коэффициентов теплопроводнос- 
ти, и на границах источников тепловыделения. Для 
погрешности менее 0,5 К и ошибки по балансу 
энергии не более 0,02 % достаточно шага по време- 
ни 0,8 мс и шага по координате 32 мкм. Получен- 
ный результат создает объективные предпосылки к 
исследованию пространственных температурных 
полей в полупроводниковых приборах. 
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